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Abstract 

The axial charge of the nucleon gA and the pion decay constant fj^ are computed in two-flavor 
lattice QCD. The simulations are carried out on lattices of various volumes and lattice spacings. 
Results are reported for pion masses as low as = 130 MeV. The volume dependence of and 
/tt can be understood quantitatively in terms of lattice ChPT. At the physical pion mass we find 
gA = 1.24(4) and = 89 it 1.1 it 1.8 MeV, using ro = 0.50(1) fm to set the scale, in good agreement 
with experiment. As a by-product we obtain the low-energy constant I4 = 4.2(1). 

PACS numbers: 12.38.Gc 
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I. INTRODUCTION 



The axial charge gA of the nucleon is a fundamental measure of nucleon structure. While 
Qa has been known accurately for many years from neutron /3 decays, a calculation of 04 
from first principles still presents a significant challenge. Present lattice calculations [l-5|, 
except perhaps underestimate the experimental value by a large amount. The resolution 
of this problem is of great importance to any further calculation of hadron structure. 

Lattice calculations of qa are in many ways connected to calculations of the pion decay 
constant f^^. Both quantities involve the axial vector current, which is not conserved and 
thus needs to be renormalized. Though it is standard practice nowadays to compute the 
renormalization constant nonperturbatively , some scope of uncertainty remains [8|. 
Another common feature is that qa and fj^ seem to be affected by large finite size corrections, 
in particular at small pion masses, which to leading order ChEFT and ChPT jo-lll appear 
to be the same in both cases. This led us to suggest to determine qa from the ratio qa/ f-K- 



First results 



12] looked indeed encouraging. 



In this paper we present our results on qa and for two flavors of nonperturbatively 
0(a) improved Wilson fermions and Wilson plaquette action |l3|. This includes simulations 
about the physical pion mass and on various lattice volumes and at various lattice spacings. 
The main focus is on finite size corrections, and the extrapolation of qa and f^^ to the 
thermodynamic limit at the physical point. 



II. LATTICE SIMULATION 



Our lattice ensembles are listed in Table HI The pion masses and the chirally extrapolated 
values of ro/a are taken from our preceding paper 13| on the nucleon mass and sigma term. 
The Sommer parameter was found to be tq = 0.50(1) fm, which we will use to set the scale 
throughout this paper. The ensembles cover three /3 values, /3 = 5.25, 5.29 and 5.40, with 
lattice spacings a = 0.076, 0.071 and 0.060 fm. 

We employ the improved axial vector current 



A^{x) = q{x)-ff,'j5q{x) + acAdf,q{x)-f5q{x) 



(1) 



where ca is taken from jl^. The improvement term does not contribute to forward matrix 



elements, but it will contribute to f-„. The calculation of qa follows 



[10, 
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with one 



/3 


K 


Volume 




9a 




ro/a 


5.25 


0.13460 


16^ X 32 


0.4932(10) 


1.442(13) 


0.0886(8) 




5.25 


0.13520 


16^ X 32 


0.3821(13) 


1.438(20) 


0.0756(8) 




5.25 


0.13575 


24^ X 48 


0.2556(5) 


1.456(10) 


0.0635(5) 


6.603(53) 


5.25 


0.13600 


24^ X 48 


0.1840(7) 


1.412(18) 


0.0550(4) 




5.25 


0.13620 


323 X 64 


0.0997(11) 


1.368(51) 


0.0439(6) 




5.29 


0.13400 


16^ X 32 


0.5767(11) 


1.437(12) 


0.0936(9) 




5.29 


0.13500 


16^ X 32 


0.4206(9) 


1.409(12) 


0.0778(5) 




5.29 


0.13550 


12^ X 32 


0.3605(32) 


1.181(60) 


0.0568(8) 




5.29 


0.13550 


16^ X 32 


0.3325(14) 


1.371(20) 


0.0675(6) 




5.29 


0.13550 


24^ X 48 


0.3270(6) 


1.459(11) 


0.0689(7) 




5.29 


0.13590 


12^ X 32 


0.3369(62) 


0.967(105) 


0.0345(9) 




5.29 
5.29 


0.13590 
0.13590 


16^ X 32 
24^ X 48 


0.2518(15) 
0.2395(5) 


1.271(32) 
1.426(7) 


0.0559(5) 
0.0588(3) 


7.004(54) 


5.29 


0.13620 


24^ X 48 


0.1552(6) 


1.334(18) 


0.0478(3) 




5.29 


0.13632 


24^ X 48 


0.1112(9) 


1.271(67) 


0.0398(4) 




5.29 


0.13632 


32^ X 64 


0.1070(5) 


1.409(24) 


0.0440(3) 




5.29 


0.13632 


40^ X 64 


0.1050(3) 


1.439(17) 


0.0445(3) 




5.29 


0.13640 


40^ X 64 


0.0660(8) 


1.363(105) 


0.0375(5) 




5.29 


0.13640 


48^ X 64 


0.0570(7) 


1.572(52) 


0.0408(11) 




5.40 


0.13500 


24^ X 48 


0.4030(4) 


1.474(7) 


0.0691(5) 




5.40 


0.13560 


24^ X 48 


0.3123(7) 


1.451(11) 


0.0620(5) 




5.40 


0.13610 


24^ X 48 


0.2208(7) 


1.410(20) 


0.0513(4) 




5.40 


0.13625 


24^ X 48 


0.1902(6) 


1.377(20) 


0.0470(3) 


8.285(74) 


5.40 


0.13640 


24^ X 48 


0.1538(10) 


1.261(34) 


0.0419(4) 


5.40 


0.13640 


32^ X 64 


0.1505(5) 


1.402(17) 


0.0442(4) 




5.40 


0.13660 


32^ X 64 


0.0845(6) 


1.206(79) 


0.0342(4) 




5.40 


0.13660 


48^ X 64 


0.0797(3) 


1.403(29) 


0.0362(3) 





TABLE I: Parameters of our lattice data sets, together with the pion mass, the bare axial charge 
and the pion decay constant. Also listed are the chirally extrapolated values of ro/a. 
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FIG. 1: The ratio R as a function of the source-sink time separation t on the 24^ x 48 lattice at 
/3 = 5.29, K = 0.13590. 

exception, namely that on the 48^ x 64 lattice at /3 = 5.29, n = 0.13640 we have employed 
Wuppertal smearing instead of Jacobi smearing. It involves computing the ratio of two- and 
three-point functions 



N and being the nucleon creation and annihilation operators at zero momentum. Any 
smearing of the source (at time 0) and sink operators (at time t) is cancelled in this ratio. 
For /3 = 5.4 we use t = 17, while the lightest two ensembles at /3 = 5.29, n = 0.13640 and 
K = 0.13632, use t = 15. All other ensembles use t = 13. 

In a ,t has bee. a.,ued that eo.tnbut.o,. f.o. exeited state, .ight be the .eason for 
qa falling short of the experimental value. To investigate this scenario, we have performed 
additional simulations on the 24^ x 48 lattice at /3 = 5.29, k = 0.13590 with a large range of 
different source-sink separations, t = 11, ■ ■ ■ , 19, albeit with somewhat lower statistics than 
our reference point at t = 13. In Fig. [1] we show the ratio R for various time separations 
t between source and sink. If true, we should find a larger value at separations t > 13. 
However, we do not see any systematic deviation of R from our result at t = 13 within the 
error bars, not even for t = 11. This provides us with confidence that our choices of t are 
sufficient with our choice of source and sink smearing. Similar conclusions were found in . 

Our smearing parameters are tuned to give a rms radius of ~ 0.5 fm, which is about half 
the radius of the nucleon. For this level of smearing no further improvement of the extracted 




(2) 
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result for qa was found by employing variational techniques [17[, which systematically sep- 
arate excited states out from the ground state at source and sink. 



The calculation of /^r follows [18[. We use the notation employed in ChPT, with the 
experimental value f.,^+ = 92.2 MeV. Our final results for the bare quantities, qa and a/^r, 
on all of our ensembles are given in Table [B 



III. FINITE SIZE CORRECTIONS 



Let us first consider the finite size corrections to qa- Utilizing the (nonrelativistic) small 
scale expansion (SSE) of the ChEFT, including pion, nucleon (N) and Z\(1232) degrees of 



freedom, we obtain to 0(e^ 



io| 



9a{oo) 



1^2 f 2 



47r2/o' 



Ki{X\n\) 
A|n| 



A{L) 



(3) 



with 



A{L) 



+ 



E 

n 

Sln'gAf^ JO 



Ko{X\n\) 



KiiMn\) 
Alnl 



dyy 



n 



KoiKy)\n\) 



Ky)\n\ 



Kii\iy)\n\) 



dyy 



/io(A(y)|n|)-^^i^i(A(y)|n|) 



(4) 



dy 



n \ / 

\n.\^0 



KiiKy)\n\) 

m\n\ 



277r/oM, 



- y 



A|n| 



where A = m^L and \{y) = f{m^,y)L with f{mT,,y) = \/rn^^ + 1/^ + 2?/Z\o, being the 
A — N mass difference. Kq and K\ denote the modified Bessel functions, and ca and g\ are 
the leading axial AN and AA couplings. The parameter ca should not be confused with 
the improvement coefficient ca in eq. ([T]). 

The second term in eq. Q, A{L\ receives contributions from chiral loops, which renor- 
malize the axial charge and act on intermediate A baryons [l^ . It turns out that the various 
contributions to A{V) effectively cancel each other over a wide range of A values. This has 
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been noticed by the authors of [19|] as welL To state an example, let us consider the 48 x 64 
lattice at /3 = 5.29, k, = 0.13640. This lattice has the lowest pion mass and is especially 
important for our final conclusions. Taking ca = 1.5 from [2^ and gi = 2.16 from SU{6), 
we find —0.044 for the total contribution, but only +0.001 for ^{L). We thus may assume 



9a{L) - gAjoo) 
9Aioo) 



Ki{Mn\) 
Alnl 



(5) 



The finite size corrections to have been computed in 111] within the context of ChPT. 
To NLO (oc m^) the outcome is 

.2 



ML) - /7r(oo) 

/7r(00) 



1^2 f 2 



47r2/, 



^i(A|n|) 
Alnl 



(6) 



The NNLO corrections are found to be very small and, thus, can safely be neglected. 

This finally shows that the leading finite size corrections to qa and f^^ are basically the 
same and will cancel in the ratio qa/ fn- Once /o, the pion decay constant in the chiral limit, 
has been fixed, expressions ([5]) and ([6]) have only one free parameter, gA{oo) and /7r(oo), 
respectively. 



The NLO correction to the pion mass reads 11 ] 



00 



00 



fi^2f2 



167r2/n 

|n|^0 



A-i(A|n|) 
A|n| 



(7) 



At smaller values of mT^L, m.,^L < 3, this expression alone cannot describe the observed 
finite size effects jisl- That is not surprising, since in a finite spatial box chiral symmetry 
does not break down spontaneously. This is because giving the system enough time it will 
rotate through all vacua. This results in a mass gap at vanishing quark masses 2ll-l23l|. 

3 



(8) 



with 



A 



+ 



1 

'0 



0.2257849591 



0.088431628 



0.8375369106 /I 



37r2 



Jin {AlL')+\n{AlL 



(9) 



where Ai are the intrinsic scale parameters of the low-ene rgy constants U = 
In (^^/"n^^phys) 1, with m^phys being the physical pion mass. In 25| we found that the 
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FIG. 2: The pion mass arriT^ as a function of lattice size for two ensembles at /3 = 5.29. The solid 
line shows a fit of eq. (fTOl) to the data. The dashed line shows the NLO result, eq. ([7]), fitted to 
the smallest mass point. 

pion mass extrapolates indeed to a finite value in the chiral limit, in good agreement with the 
expected resul t (El) . This also has an effect on m^r in the region of small, but nonvanishing, 
25|. We thus expect the finite size correction to be effectively given by 



quark masses 



/T\ / N \ ^ A"i(A|n|) 3c(m7r) 



2 f2 
JO n 
|n|^0 



with the parameter c{mT^) rapidly dropping to zero at larger pion masses. 



(10) 



IV. EXTRAPOLATION TO INFINITE VOLUME 



In the following fits we take /o = 86 MeV 26[ . There is some freedom in the value of 
the pion mass to take in eqs. (E]), (E]) and ffTU]) . We choose = 171^,(00) in A, X{y) and 
c(m^), and m^r = mT^{L) otherwise. 

Let us first consider the pion mass. In Fig. |2] we show the fits of eq. (fTO]) to for two 
of our lattice ensembles. The corrections to m,r are well described by this equation. Apart 
from m7r(oo), we have one free parameter, c(m^), only. Equally good fits are obtained for 
/3 = 5.40, K = 0.13660 and 0.13640. The parameter c{m.^) is found to vanish with a large 
inverse power of the pion mass. The finite size corrections predicted by the NLO expression 
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FIG. 3: The bare axial charge qa and the bare pion decay constant 0/,^ as a function of the spatial 
extent of the lattice, together with the leading order finite size corrections of eqs. dSJ and ([B]). 
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([7]), on the other hand, are nowhere near as big as the effect shown by the data. In Table ITTl 
we hst our final pion masses. Our lowest mass turns out to be = 130(5) MeV. 

Let us now turn to the axial charge and the pion decay constant. In Fig. |3]we show the 
fits of eqs. and (jH]) to qa and a/^, respectively, for our three lowest pion masses. It shows 
that the finite size effects are well described by the leading order formulae. The corrections 
are large, in particular at the physical pion mass, displayed by the top two figures. On the 
L = 3 fm lattice the corrections are of the order of 10% still. 

To obtain continuum numbers, we need to renormalize the axial vector current. That 
reads 

= ZA{l + bAamg) A^. (11) 

The coefficient Ba is required to maintain 0{a) improvement for nonvanishing quark masses 
niq as well. The renormalization constant Za has been computed nonperturbatively in [3] , 
employing the Rome-Southampton method 6|, with the result 

(3 5.25 5.29 5.40 , ^ 

(12) 

Za 0.760(1) 0.764(1) 0.777(1) 

The coefficient Ba is only known perturbatively j^^], i.e. 

5^ = 1 + 0.1522/. (13) 

We denote the renormalized axial charge and pion decay constant by Qa and f^, respectively. 

In Table Hn we give our extrapolated numbers. For pion masses < 300 MeV we demand 
that we have at least two lattice volumes to ensure a controlled extrapolation. This excludes 
the point at /3 = 5.25, k = 0.13620 with pion mass ^ 250 MeV. 



V. FINAL RESULTS AND DISCUSSION 

Our results for are plotted in Fig. |H Since after finite volume corrections our lightest 
pion mass is 130 MeV, no extrapolation is required. Instead, we quote our result at our 
lightest pion mass, namely 

= 1.24 ±0.04, (14) 

in good agreement with the experimental value. It turns out that Qa hovers around ^1.1 
for > 250 MeV, a feature it shares with most other lattice calculations Only on the 

10 



p 


K 


[MeV] 


9a 


[MeV] 


5.25 


0.13600 


479(2) 


1.07(1) 


108.9(0.8) 


5.29 


0.13620 


426(2) 


1.05(2) 


103.6(0.6) 


5.29 


0.13632 


284(2) 


1.10(2) 


94.7(0.6) 


5.29 


0.13640 


130(5) 


1.24(4) 


89.7(1.5) 


5.40 


0.13640 


492(2) 


1.09(1) 


112.3(0.9) 


5.40 


0.13660 


253(2) 


1.09(2) 


93.0(0.7) 



TABLE II: The pion mass and the renormahzed axial couphng and pion decay constant extrapo- 
lated to the infinite volume for m-,^ < 500 MeV and ro = 0.50 fm. 



very last 100 MeV from the physical point does rise to its final value. 



is not totally unexpected, from general arguments 



28 



and from ChEFT 



his phenomenon 
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the chiral limit ChEFT predicts, following the notation of 



io| 



30|. Near 



gV2 



ml + 4: [-B9(m^phys) - 2^1^ 52o(m^phys 



(15) 



ml ln(m^/m^phys) + O (ml) 
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29] and 30|], are equivalent with Bg = die 



To this order, both sets of chiral expansions, 
and B20 = d28- In (fTSjl we have chosen A = m,rphys (A being the scale parameter of the 

dimensional regularization) . The coupling i?2o cannot be observed independently of Bg. 

-2 



Taking -820(^^71 phys) = 0, the preferred value of B^ is [311] -Bg ('t^tt phys 



-1.4 ±1.2) GeV~ 



A fit of the leading order chiral formula 0151) to the data points in Fig. H] is shown by the 
shaded area. The fit gives g'^ = 1.26(7) and -B9(m^phys) = (—2.1 ± 1.0) GeV~^. We do not 
see any scaling violations in the region where we have results at multiple lattice spacings, 
i.e. m^r > 250 MeV. However, below that region we cannot draw any further conclusions. 

Our results for are plotted in Fig. [5l Again, no extrapolation to the physical point is 
needed, where we obtain from our lightest simulation point 



= 89.7 ± 1.5 ± 1.8 MeV; 



(16) 



using To = 0.50(1) fm. The second error in eq. f|T6|) is due to the error on tq. In the case of 
we do not see any scaling violations either. 
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FIG. 4: The renormalized axial charge in the infinite volume plotted against m^(oo), together 
with the experimental value qa = 1-27 (x). The shaded area shows the fit of eq. (jlSp to the data. 
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FIG. 5: The renormalized pion decay constant in the infinite volume plotted against m^(oo), 
together with the experimental value fj^ = 92.2 MeV (x). The curve shows a fit of eq. (jl7p to the 
data. 
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Instead of taking from the entry at the lowest pion mass, eq. (fT6|) . it might be a better 
idea to include the adjacent data points in the analysis as well and fit the data by a chiral 
ansatz |2J], 



fn — fo 



2 

ml 



In (Al/ml) 



Ami. (17) 



TV 

The result of the fit is shown in Fig. [51 At the physical point we obtain the result 

= 89.6 ± 1.1 ± 1.8 MeV, (18) 

which is fully consistent with the number in eq. f lT6|) . The main effect is that the statistical 
error has reduced by 30%. In the chiral limit we obtain /q = 86(1) MeV, which agrees with 
the number used in the finite size correction formulae. A fit of the chiral ansatz ( IT71) to the 
lowest four data points with A = gives the low-energy constant 

h = ln{Al/ml^^^^)=A.2±0A. (19) 

Though eqs. (|5]) and (|6]) describe the finite size corrections very well, it should be noted 
that other choices of the pion decay constant (for example instead of /o) are possible. 
There is also some ambiguity in the determination of the renormalization constant Za js]. 
Finite size corrections and uncertainties of Za can be evaded by considering the ratio qa/ fn- 
In Fig. [6] we plot the ratio Qa/ fn for our (raw) data points listed in Table HI restricting 
ourselves to pion masses < 750 MeV and again taking tq = 0.50(1) fm to set the scale. If 
we have more than one volume at a given n value, we show the result of the largest volume. 
The data points for all three /3 values lie nicely on a universal curve. The lowest pion mass 
in Fig. [6] is 157 MeV (as opposed to 130 MeV in the infinite volume). 

The leading order chiral expansion of Qa/ fn (in the infinite volume) can be cast in the 
form Q, y, Q 

^-^ = A + Bml + Cml\nml + Dml. (20) 

We have fitted eq. (l20l) to the data points in Fig. [5l The result is shown by the solid curve. 
At the physical point this gives 

= 13.95 ± 0.71 ± 0.30 GeV-^ (21) 

Jtt 

Again, the second error is due to the error on tq. Multiplying the ratio (12T]) by the physical 
value of /tt, = 92.2 MeV, we then obtain 

= 1.29 ±0.05 ±0.03. (22) 

13 




FIG. 6: The ratio qa/ f-K-, together with the experimental value (x). The curve shows a fit of 
eq. ([201) to the data. 

This result is in agreement with the direct calculation of g^^ given in eq. (fT4l) . which involved 
extrapolation of qa to the infinite volume and nonperturbative renormalization. 

In Fig. iniwe have set the scale by Tq = 0.50(1) fm. Alternatively, we could have set the 
scale by the physical value of f^^. That would give the value 

= 1.27 ±0.05. (23) 
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